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B-hadrons are an ideal tool for advancing our current understanding of the flavor sector of the Standard Model 
(SM). The study of B-meson production and decays is one of the key physics themes at the Large Hadron 
Collider (LHC) thanks to the large production rate and the fact that B-hadrons are relatively easy to trigger on 
and identify due to their long lifetime and decays to muons. This talk presents the cross section measurement 
for the exclusive final state B^ — >■ J/if> (f> and an evaluation of the decay branching fraction from the previously 
published exclusive-B production cross sections. Besides probing the heavy quark properties for the first time 
at the LHC energy, these measurements are also important tools for understanding and calibrating the detector, 
giving input for Monte Carlo tuning, and providing results for direct comparison with other experiments. 

o 

^ 1. Heavy-flavor physics at the CMS experiment 

C/3 : 

On : The measurements of differential cross sections for heavy-quark production in high-energy hadronic inter- 
actions are critical input for the underlying next-to-leading order (NLO) Quantum Chromodynamics (QCD) 
,— i . calculations Q| . While progress has been achieved in the understanding of heavy-quark production at Tevatron 
>^ ■ energies 10j, large theoretical uncertainties remain due to the dependence on the renormalization and fac- 
torization scales. Measurements of b-hadron production at the higher energies provided by the LHC represent 
O ■ an important new test of theoretical approaches that aim to reduce the scale dependence of NLO QCD calcula- 
tions Jll|, [HJ]. The Compact Muon Solenoid (CMS) experiment, that covers a rapidity range complementary to 
i-^ ■ the specialised b-physics experiment LHCb [13|, recently measured the cross sections for production of B + 14 1 
and B°[l^] in pp collisions at y/s = 7 TeV. This talk presents the first measurement of the production of B", 
H ■ with B° decaying into J/i/j <fi, that adds information to the understanding of b-quark production at this energy. 
' Data and theoretical predictions are compared to NLO predictions of heavy-quark production. The result, 
combined with the published B + and B° production cross section measurements, is then used for the evaluation 
**** ■ of the B° — > J/ip 4> branching fraction. The cross section measurement for the B^ — > J ftp (f> decay is only the 



o 



- 1—1 

X 



first part of a longer term program that includes the measurement of the B s mesons properties (lifetime and 
lifetime difference) and CP-violating parameters (CP-even and CP-odd amplitude strenghts and weak phase). 



2. The CMS detector 



A detailed description of the CMS detector can be found elsewhere [16J. The primary components used in 
this analysis are the silicon tracker and the muon systems. The tracker operates in a 3.8 T axial magnetic field 
generated by a superconducting solenoid having an internal diameter of 6 m. The tracker consists of three 
cylindrical layers of silicon pixel detectors complemented by two disks in the forward and backward directions. 
The radial region between 20 and 116 cm is occupied by several layers of silicon strip detectors in barrel and 
disk configurations, ensuring at least nine hits in the pseudorapidity range |?7| < 2.4, where rj = — In [tan {0/2)] 
and is the polar angle of the track measured from the positive z-axis of a right-handed coordinate system, with 
the origin at the nominal interaction point, the a;- axis pointing to the centre of the LHC, the y-axis pointing 
up (perpendicular to the LHC plane), and the z-axis along the counterclockwise-beam direction. An impact 
parameter resolution of about 15 yum and a pt resolution of about 1.5% are achieved for charged particles 
with transverse momenta up to 100 GeV/c. Muons are identified in the range |?7| < 2.4, with detection planes 
made of drift tubes, cathode strip chambers, and resistive plate chambers, embedded in the steel return yoke. 
The first level of the CMS trigger system uses information from the crystal electromagnetic calorimeter, the 
brass/scintillator hadron calorimeter, and the muon detectors to select the most interesting events in less than 
1 /is. The high level trigger (HLT) employs software algorithms and a farm of commercial processors to further 
decrease the event rate using information from all detector subsystems. The events used in the measurement 
reported in this paper were collected with a trigger requiring the presence of two muons at the HLT, with no 
explicit momentum threshold. 
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3. Muon reconstruction 

In CMS muons are denned as tracks reconstructed in the silicon trackers and associated to a compatible 
signal in the muon chambers. Two different muon types are available in CMS. The first one, referred to as a 
Global Muon, provides high-purity reconstruction for muons with px > 4 GeV/c in the central pseudo-rapidity 
region |nj < 1.5, and px > 1 GeV/c in the forward region. Global Muons are built as a combined fit of silicon 
and muon-chamber hits, belonging to independent tracks found in the tracker and muon systems. The second 
muon type, referred to as a Tracker Muon, achieves a better reconstruction efficiency at lower momenta. The 
requirements for a Tracker Muon are relaxed compared to the Global Muons, at the expense of a slightly larger 
background: tracks found in the tracker matched to only one muon segment are accepted and not refitted. If 
two (or more) tracks are close to each other, it is possible that the same muon segment or set of segments is 
associated to more than one track. In this case the best track is selected based on the matching between the 
extrapolated track and the segment in the muon detectors. 



4. Strategy outline 

A sample of exclusive B° — > J/ip 4> decays, with J/ip —> /i + p~ and 4> - > K + K~ , is reconstructed from the data 
collected in 2010 by the CMS experiment, corresponding to an integrated luminosity of 39.6 ± 1.6 pb -1 . The 
differential production cross sections, da/dp B and da/dy B , are determined as functions of the transverse mo- 
mentum p B and rapidity |y B | of the reconstructed B° candidate. Here, the rapidity y B is defined as \ In |jrf^S 
where E is the particle's energy and pl is the particle's momentum along the counterclockwise beam direction. 
The differential cross sections are calculated from the measured signal yields (n s j g ), corrected for the overall 
efficiency (e), bin size (Ax, with x = p B , |y B |), and integrated luminosity (L) as 

Ar(pp -» B° -s- J/-0 (f>) = n sig 

dx 2-e-B-L-Ax ' [ ' 

where B is the product of the branching fractions for the decays of the J/ip and <j> mesons. In each bin the 
signal yield is extracted with an unbinned maximum likelihood fit to the J ftp (f> invariant mass and proper decay 
length ct of the B° candidates. The factor of two in Eq.[T]is required since we report the result as a cross section 
for Bg production alone, while both B° and B^ are included in n s - lg . The size of the bins is chosen such that the 
statistical uncertainty on n s i g is comparable in each of them. 



4.1. Reconstruction and B° selection 

Reconstruction of B° — > J/ip 4> candidates begins by identifying J/ip — > decays. The muon candidates 

must have one or more reconstructed segments in the muon system that match the extrapolated position of a 
track reconstructed in the tracker. Furthermore, the muons are required to lie within a kinematic acceptance 
region defined as: p£ > 3.3 GeV/c for |<| < 1.3; total momentum p^ > 2.9 GeV/c for 1.3 < |r? M | < 2.2; and 
> 0.8 GeV/c for 2.2 < < 2.4. Two oppositely charged muon candidates are paired and are required to 
originate from a common vertex using a Kalman vertex fit [171 ] . The muon pair is required to have a transverse 
momentum px > 0.5 GeV/c and an invariant mass within 150 MeV/c 2 of the world average i /ip mass value [l8| . 
which corresponds to more than three times the measured dimuon invariant mass resolution [19(. 

Candidate <f> mesons are reconstructed from pairs of oppositely charged tracks with px > 0.7 GeV/c that are 
selected from a sample with the muon candidate tracks removed. The tracks are required to have at least five 
hits in the silicon tracker detectors, and a track \ 2 per degree of freedom less than five. Each track is assumed 
to be a kaon and the invariant mass of a track pair has to be within 10 MeV/c 2 of the world average 0-meson 
mass fl8l |. 

The B^ candidates are formed by combining a J /tp candidate with a <j> candidate. The two muons and the two 
kaons are subjected to a combined vertex and kinematic fit [2p| . where in addition the dimuon invariant mass 
is constrained to the nominal J/ip mass. The selected candidates must have a resulting \ 2 vertex probability 
greater than 2%, an invariant mass between 5.20 and 5.65 GeV/c 2 , and be in the kinematic range 8 < p B < 
50 GeV/c and \y B \ < 2.4. For events with more than one candidate, the one with the highest vertex-fit 
probability is selected. This results in the correct choice 97% of all cases, as determined from simulated signal 
events. 

The proper decay length of each selected B^ candidate is calculated using the formula ct = c(Mb/p b )L xv , 
where the transverse decay length L xy is the length of the vector s pointing from the primary vertex [2l| to the 
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secondary vertex projected onto the Bj? transverse momentum: L xy — (s ■ p^)/p®, with Mb the reconstructed 
mass of the Bg candidate. Candidate B° mesons are accepted within the range —0.05 < ct < 0.35 cm. A total 
of 6, 200 events pass all the selection criteria. 

The efficiency of the B^ reconstruction is computed with a combination of techniques using the data and 
large samples of simulated signal events generated using PYTHIA 6.422 [22J. The decays of unstable particles 
are described by the EVTGEN [23| simulation. Long-lived particles are then propagated through a detailed 
description of the CMS detector based on the GEANT4 [24J package. The trigger and muon-reconstruction 
efficiencies are obtained from a large sample of inclusive J/ip —> decays. The total efficiency of this 

selection, defined as the fraction of B[? — > J/-0 <fi decays produced with 8 < < 50 GeV/c and |i/ B | < 2.4 that 
pass all criteria, ranges from 1.3% for p^ w 8 GcV/c to 19.6% for p® > 23 GeV/c. 



4.2. Muon efficiency 

The "Tag and Probe" [l^| technique is a data driven method to measure the single muon tracking, identifi- 
cation and trigger efficiencies. It makes use of a well-known dimuon resonance (such as J/i/> mesons) to supply 
tags and probes. The choice of such resonance is due to the CMS experiment ability to measure muon momenta 
with high precision and reconstruct and identify muons with high efficiencies. Events are selected with strict 
selection requirements on one muon (tag), and with a more relaxed selection on the other track (probe), such 
that the selection applied to the probe track does not bias the efficiency that one wants to measure. The probe 
tracks are separated into two categories depending on whether they pass or fail the more restrictive selection. 
If Ptag and Pprobe are the four-momenta of the tag muon and probe track respectively then the invariant mass 
m of the particle is given by: 

m = yjptag + Pprobe (2) 

The technical chain runs as follows: all events of the J/ip samples are passed through. An event is kept if a 
tag muon + probe track combination is found while satisfying the predefined criteria. In this case the invariant 
mass of the combination enters a muon-track mass histogram. The criteria for the choice of both, tag and probe, 
should be defined such that background is reduced and the J/tp mass peak is clearly visible. If the track is 
subsequently identified as a global/tracker muon the invariant mass of the tag and probe combination enters 
another muon-muon mass histogram. The latter should contain almost exclusively J/ip — > events. A 

simultaneous unbinned Maximum Likelihood (ML) fit to both mass distribution extracts the efficiency. The 
dimuon efficiencies are calculated as the product of the single-muon efficiencies obtained with this method. 
Corrections to account for correlations between the two muons (1-3%) are obtained from simulation studies. 
The correction factors are determined in bins of single muon p!^ and rf 1 , and are applied independently to each 
muon from a B° — > 3 ftp <f> decay in the simulation to determine the total corrected efficiency. 



5. Fit technique 

The two main background sources are prompt and non-prompt J/ijj production. The latter background is 
mainly composed of B + and B° mesons that decay to a 3/tp and a higher-mass K-meson state (such as the ). 
Such events tend to have lower reconstructed Mb mass. Inspection of the reconstructed J/ip4> invariant mass 
for a large variety of potential B background channels confirms that there is no single dominant component 
and that the channel B° J/ip K*° (with K*° K + n~), which a priori is kinematically similar to the signal 
decay and more abundantly produced, is strongly suppressed by the restriction on the K + K~ invariant mass. 
A study of the sidebands of the dimuon invariant mass distribution confirms that the contamination from events 
without a J/ip decay to two muons is negligible after all selection criteria have been applied. 

The signal yields in each p^ and |y B | bin are obtained using an unbinned extended maximum-likelihood fit to 
Mb and ct. The likelihood for event j is obtained by summing the product of the yield ni and the probability 
density functions (PDF) Vi and Qi for each of the signal and background hypotheses i. Three individual 
components are considered: signal, non-prompt b — > 3/?p X, and prompt J/ip. The extended likelihood function 
is then the product of likelihoods for each event j: 



C = exp 



■5> n 



i=l 



riiPi(M B ; ai)Qi(ct; A) 



(3) 
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The PDFs Vi and Qi are parameterized separately for each fit component with shape parameters c?j for Mb 
and Pi for ct. The yields ni are then determined by minimizing the quantity — ln£ with respect to the signal 
yields and a subset of the PDF parameters [lf|. The PDFs are constructed from basic analytical functions 
that satisfactorily describe the variable distributions from simulated events. Shape parameters are obtained 
from data when possible. The Mb PDF is the sum of two Gaussian functions for the signal, a second-order 
polynomial for the non-prompt J/tp that allows for possible curvature in the shape, and a first-order polynomial 
for prompt J/i/j. The resolution on Mb is approximately 20 MeV/c 2 near the Bj? mass. 

For the signal, the ct PDF is a single exponential parameterized in terms of a proper decay length ct. It is 
convolved with a resolution function that is a combination of two Gaussian functions to account for a dominant 
core and small outlier distribution; the core fraction is varied in the fit and found to be consistently larger than 
95%. The ct distribution for the non-prompt J/ip background is described by a sum of two exponentials, with 
effective lifetimes that are allowed to be different. The "long-lifetime exponential" corresponds to decays of b- 
hadrons to a 3/tp plus some charged particles that survive the cf> selection, while the "short-lifetime exponential" 
accounts for events where the muons from the J /ip decay are wrongly combined with hadron tracks originating 
from the pp collision point. The exponential functions are convolved with a resolution function with the same 
parameters as the signal. For the prompt J/V> component the pure resolution function is used. The core 
resolution in ct is measured in data to be 45 /xm. All background shapes are obtained directly from data, while 
the signal shape in Mb is taken from a fit to reconstructed signal events from the simulation. The effective 
lifetime and resolution function parameters for prompt and non-prompt backgrounds are extracted from data in 
regions of Mb that are separated by more than four times the width of the observed B° signal from the mean B° 
peak position (M B sidebands): 5.20 < M B < 5.29 GeV/c 2 and 5.45 < M B < 5.65 GeV/c 2 . A comparison of 
the PDF shapes for the different sideband regions in simulated events confirms that they are extrapolated well 
into the signal region. With the lifetimes for signal and non-prompt background fixed from this first step, the 
resolution function parameters are then determined separately in each p B and |y B | bin, from the Mb sidebands. 
The signal and background yields in each p B and \y B \ bin are determined in a final iteration, using the full Mb 
range, with all parameters floating except the background lifetimes and the lifetime resolution functions, which 
are fixed to the results of the fit to the Mb sidebands. It has been verified that leaving all parameters floating 
changes the signal yield by an amount smaller than the systematic uncertainty assigned to the fit procedure. 

5.1. Fit results 

Figure Q] shows the fit projections for Mb and ct from the inclusive sample with 8 < p® < 50 GeV/c and 
|y B | < 2.4. When plotting Mb, the selection ct > 0.01 cm is applied for better visibility of the individual 
contributions. The number of signal events in the entire data sample is 549 ± 32, where the uncertainty 
is statistical only. The obtained proper decay length of the signal, ct = 478 ± 26 /im, is within 1.4 standard 
deviations of the world average value jig , even though this analysis was not optimized for lifetime measurements. 



6. Systematics 

The cross section measurement is affected by several sources of systematic uncertainty arising from uncertain- 
ties on the fit, efficiencies, branching fractions, and integrated luminosity. In every bin the total uncertainty is 
about 11%. Uncertainties on the muon efficiencies from the trigger, identification, and tracking are determined 
directly from data (3 — 5%). The uncertainty of the method employed to measure the efficiency in the data has 
been estimated from a large sample of full-detector simulated events (1 — 3%). The tracking efficiency for the 
charged kaons is consistent with simulation. A conservative uncertainty of at most 9% in each bin has been 
assigned for the hadronic track reconstruction (adding linearly the uncertainties on the two kaon tracks (25j). 
which includes the uncertainty due to misalignment of the silicon detectors. The uncertainty of the fit pro- 
cedure arising from potential biases and imperfect knowledge of the PDF parameters is estimated by varying 
the parameters by one standard deviation (2 — 4%). The contribution related to the B^ momentum spectrum 
(I — 3%) is evaluated by reweighting the shape of the p B distribution generated with PYTHIA to match the 
spectrum predicted by MC@NLO [27|. An uncertainty of 1% is assigned to the variation of the selection criteria 
applied to the vertex fit probability, the transverse momentum of the kaons, the B^ transverse momentum, and 
the K + K~ invariant mass. An uncertainty is added to account for the limited number of simulated events 
(at most 3% in the highest p B bin). The total uncorrelated systematic uncertainty in each bin is the sum in 
quadrature of the individual uncertainties (see Table HJ . In addition, there are common uncertainties of 4% 
from the integrated luminosity measurement [28| and 1.4% from the 3/tp and <fi branching fractions 18]. As 
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Figure 1: Projections of the fit results in Mb (a) and ct (b) for 8 < < 50 GeV/c and \y B \ < 2.4. The curves in 
each plot are: the sum of all contributions (solid line); signal (dashed); prompt 3/t[> (dotted); and non-prompt J ftp 
(dot-dashed). For better visibility of the individual contributions, plot (a) includes the requirement ct > 0.01 cm. 



the reported result is a measurement of the Bj? cross section times the — > i/ip <fi branching fraction, the 30% 
uncertainty in the B° — > J/ip <f> branching fraction [18( is not included in the result. 
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Table I: Summary table of the relative systematic uncertainties in the measurement of the B° — > J/i/> <f) production cross 
section. 



Source 


TJnccrtciiiity (%) 


Ivluon Reconstruction Efficiency 


o — 


1— T q r\ r r\n ' 1 V - ufl/in rr TVtti nion /"m r 
1 J.dLl.1 Ull _LI £tCft.lIlJi J2jllH_,lt;IH_, y 


7 8 
l .0 


R ppnn crt"rn o\ i nn Tltti c*\ PTifv 


2 — 3 


1VT i t; a 1 1 cm in pti f 


2 — 4 


PT~|y 1 Spectrum 


1-3 


Probability Density Function 


2-4 


Uncorrelated Systematic Errors 


10 - 11 


Branching Fractions 


1.4 


Luminosity 


4 


Correlated Systematic Errors 


4.2 


Total Systematic Error 


11 - 12 



7. Differential cross section measurement 



The differential cross sections times branching fraction as functions of p B and | y B | are plotted in Fig. [2j 
together with predictions from MC@NLO and PYTHIA . The predictions of MC@NLO use the renormalization 
and factorization scales fi — y^m^c 4 + p\c 2 , where pt is the transverse momentum of the b quark, a b-quark 
mass of mb = 4.75 GeV/c 2 , and the CTEQ6M parton distribution functions j29j. The uncertainty in the 
MC@NLO cross section is obtained simultaneously varying the renormalization and factorization scales by 
factors of two, varying nib by ±0.25 GeV/c 2 , and using the CTEQ6.6 parton distribution function set. The 
prediction of PYTHIA uses the CTEQ6L1 parton distribution functions [29| . a b-quark mass of 4.8 GeV/c 2 , 
and the Z2 tune [3(| to simulate the underlying event. The total integrated B° cross section times B° — > J/-0 <f) 
branching fraction for the range 8 < p® < 50 GeV/c and |y B | < 2.4 is measured to be (6.9 ±0.6 ±0.6) nb, where 
the first uncertainty is statistical and the second is systematic. The statistical and systematic uncertainties are 
derived from the bin-by-bin uncertainties and propagated through the sum. The measured total cross section 
lies between the theoretical predictions of MC@NLO (4.6±\-% ± 1.4 nb) and PYTHIA (9.4 ± 2.8 nb), where 
the last uncertainty is from the B° — > J/?/> <fi branching fraction [18j . Also the previous CMS cross-section 
measurements of B + 14 1 and B° [T5| production in pp collisions at y^i = 7 TeV, gave values between the two 
theory predictions, indicating internal consistency amongst the three different B-mcson results. 



8. Evaluation of the branching fraction for B° — >- J/?/> <fi 

The branching fraction BF(B° — » J/tp <j>) can be calculated independently with respect to either the B + — > 
J/-0K+ or the B° — > J/ipK° decay channels. The observed cross section for the decay mode B° — >• J/-0 cf) can 
be written as: 

a(pp -J- B° -j. J/V» 0) - o{pp >/.:../;• BF(B° S J/^ 0) • /£, (4) 
and similarly, for the B + and B° mode: 

a(pp -> B+A) = <t(pp -> 6) • /„ • (5) 
a(pp -> B°A) = a( W -> 6) • / rf • / B ° n (6) 

Here, the /„, fd, and / s are the probabilities that the b anti-quark will hadronize and form a B + , B°, and B^ 
meson, respectively. The fractions / B n correct for the limited range in rapidity and transverse momentum in 
the different analyses. The extrapolation to the full kinematic range is theory dependent. The NLO theory 
predictions [27j for the expected differential cross section [Til [l5], [3l[ values are in good agreement with the 
measured ones in each of the three decay channels. Therefore, it is possible to identify the best central model 
(CTEQ6M, Qji — Qp = 1, and mb — 4.75 GeV/c 2 ) to predict the full kinematic range in p B and |y B |. The 
spectra are obtained from a number of NLO generated events large enough to keep the relative statistical error 
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Table II: Summary table for the fragmentation fractions used in the evaluation of the B° — > J/ip 4> branching fraction. 



LEP+Tevatron 


r(6 -)• B s °) 
r(6 ->■ b+) 
r(6^B°) 


(11.0 ±1.2)% 
(40.3 ±1.1)% 
(40.3 ±1.1)% 


fs _ r(B») 
fu.d ~ r(B+-") 


(27.2 ±3.1)% 


Tevatron 


r(5 ->• b°) 
r(6 ->■ b+) 
r(6^B°) 


(11.1 ±1.4)% 
(33.9 ±3.1)% 
(33.9 ±3.1)% 


/« r(B^) 
/«,d - r(B+.») 


(32.7 ±5.1)% 



on the predicted ratio at 0.1%. The values for the branching fractions are derived from the measured cross 
sections: 



BF(B° 3/^ 0) 



a(pp 



B°->J/^) f Utd f£ 



a(pp -> B+>°X) 



J k. 



BO 
kin 



(J) 



Both values depend on the ratio of the fraction of kinematic ranges that vary slowly over a wide range of model 
parameters. The values for the fragmentation fractions are obtained from HFAG, printed in PDG [18| under b- 
hadron admixtures. Both results, from Tevatron and from the combination of LEP and Tevatron measurements, 
are considered and reported in Table HU The sources of uncertainties are listed in Table [ 



Table III: Summary table of the systematic uncertainties for the evaluation of the B° — > J /if) <f> branching fraction. 



Source 


R P B s^J/</> 
Dr B+->J/,/.K+ 






Experimental Uncertainties 


Cross Section 


15.8 


16.5 


NLO Spectrum 


4.6 


4.3 




PDG Uncertainties 


Branching Fractions 


3.5 


3.8 


Fragmentation Fractions 


11.2 


11.2 



The experimental error is the sum in quadrature of the statistical and systematic uncertainties calculated 
as described in the cross section measurement papers [l4], EH HH . It contains the uncertainty of the yield as 
extracted from maximum likelihood fits, the uncertainty on the reconstruction and hadron-tracking efficiencies, 
misalignment, and the variation with different models to correct for the limited kinematic ranges. From that, we 
distinguish the uncertainties of the branching fractions and the fragmentation fractions in PDG [l8| . Individual 
uncertainties are added in quadrature. The cross section measurements in the different exclusive B decay modes, 
omitting the uncertainty due to the luminosity measurement, are: 

• a{pp B° -> J/tp 4>) = (6.9 ± 0.6 ± 0.5) x lO" 3 ^ 

• a(pp -> B+X) = (28.3 ± 2.4 ± 2.0)^b 

• a(pp -> B°X) = (33.2 ± 2.5 ± 3.1)pb 

The branching fraction BF(Bg — > 3 /ip (p) is calculated according to Eq. [7] with respect to the B + meson and to 
the B° meson productions. The two measurements of the branching fraction BF(Bj? — » 3/ip cp) overlap within 
one standard deviation. Their error weighted average is 



BF(B° -> 3/ip (p) = (1.8 ± 0.2 ± 0.2) x 10" 3 



(8) 
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where the errors are the combined experimental and PDG uncertainties, respectively. The result is listed in 
the same way as PDG does for the present world average. The result agrees within one standard deviation 
with the PDG value [H[ of BF(B° — !• J/ip 4>) = (1.4 ± 0.4 ± 0.2) x 10~ 3 assuming the same uncertainty on the 
b-quark fragmentation fraction. Using the fragmentation fractions as extracted from the Tevatron, only, the 
error weighted average for the branching fraction is 

BF(B° -> J/Vj 4>) = (1.5 ± 0.2 ± 0.2) x 10~ 3 (9) 

to be compared with the PDG value BF(B° J/tp 4>) = (1.2 ± 0.3 ± 0.2) x 10~ 3 that has been extracted in the 
same way. 

9. Summary 

In summary, the first measurements of the B° differential cross sections da /dp® and da/dy B in pp collisions at 
^/s = 7 TeV and in the decay channel Bg — > J /if) 4> have been presented. The measurement has been performed 
in four bins in the kinematic range \y B \ < 2.4 and 8 < p® < 50 GeV/c. This study complements previous 
results in moving towards a comprehensive description of b-hadron production at \/s = 7 TeV. An estimation of 
the B° — > J/ip (j> branching fraction has been calculated from the published CMS measurements of inclusive B + 
and B° production cross section. We calculate the branching fraction for the B° — > J/i(j <f> decay, assuming the 
fragmentation fractions extracted from measurements at the Tevatron, to (1.5 ± 0.2 ± 0.2) x 10~ 3 that agrees 
within one standard deviation with the value published from the Tevatron. 
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Figure 2: Measured differential cross sections da/dp^ (a) and da/dy B (b) compared with theoretical predictions. The 
(yellow) band represents the sum in quadrature of statistical and systematic uncertainties. The dotted (red) line is the 
PYTHIA prediction; the solid and dashed (blue) lines are the MC@NLO prediction and its uncertainty, respectively. 
The common uncertainties of 4% on the data points, due to the integrated luminosity, and of 30% on the theory curves, 
due to the B° — > J/ip <f> branching fraction, are not shown. 
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